COMPLEX RANDOM MATRIX MODELS WITH 
POSSIBLE APPLICATIONS TO SPIN-IMPURITY 
SCATTERING IN QUANTUM HALL FLUIDS 



S. Hikami 1 and A. Zee 2 

Department of Pure and Applied Sciences 1 
University of Tokyo, Meguro-ku, Komaba, Tokyo 153, Japan 

Institute for Theoretical Physics 2 
University of California Santa Barbara, CA 93106, USA 

Abstract 

We study the one-point and two-point Green's functions in a com- 
plex random matrix model to sub-leading orders in the large TV" limit. 
We take this complex matrix model as a model for the two-state scat- 
tering problem, as applied to spin dependent scattering of impurities 
in quantum Hall fluids. The density of state shows a singularity at 
the band center due to reflection symmetry. We also compute the one- 
point Green's function for a generalized situation by putting random 
matrices on a lattice of arbitrary dimensions. 
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1 Introduction 



Recently, matrix modes have been studied in various contexts. In particular, 
in a series of papers [1-4], the universality of the connected two-point corre- 
lation function was discussed. These studies focussed on N by iV Hermitian 
matrices with TV approaching infinity. 

In this paper, we investigate the behavior of the Green's functions of 
a complex matrix model. To leading order in the large N limit, complex 
matrix model and Hermitian matrix model behave similarly. For instance, 
the density of state, for Gaussian randomness, obeys Wigner's semi-circle 
law in both cases. However, in subleading orders, complex random matrix 
and Hermitian random matrix behave differently. In the language of random 
surfaces, a random complex matrix model represents, due to the complex 
conjugate pairs, a surface made of plaquets with two colors, like a red and 
black checkerboard, with the rule that a given plaquette can only be glued 
to a plaquette of a different color. 

Our motivations stem from possible applications of the complex random 
matrix model to physical problems involving scattering on impurities. In 
particular, recently Hikami, Shirai, and Wegner [5] has proposed a model for 
impurity scattering in quantum Hall fluids in the spin degenerate case. For 
certain quantum Hall samples, disorder broadening can be much larger than 
the Zeeman splitting between spin up and spin down electrons, in which case 
the spin up and spin down electrons have the same energy. In [5] the further 
simplification is made that when a spin up electron scattering on an impurity 
becomes a spin down electron, and vice versa. This is known as the "strong 
spin orbit" case, in which case it is known that an extended state appears 
at the band center of the lowest Landau level with white noise Gaussian 
random scattering and that the density of state shows a singularity at the 
band center [5-9]. Although this problem has been simulated numerically [9], 
the nature of the singularity of the density of state remains unclear. Another 
example involving scattering between two states occurs in high temperature 
superconductors, in which the conducting plane contains two different sites, 
copper and oxygen. The density of state also shows a singular behavior at 
the band center. 

More generally, the problem of scattering between C sectors (with C = 2 
in the example mentioned above) is of interest. It turns out that a class of 
matrix models called "lattices of matrices" and studied by Brezin and Zee 
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[3] is relevant for this class of problem. In these models, random matrices 
are placed on a lattice of arbitrary dimension [10]. It was shown that in the 
large N limit, various correlation functions can be determined. The two- 
state scattering problem considered here corresponds to the simple case of a 
lattice consisting of two points. Following the analysis of [3] we can readily 
generalize the two- state scattering problem to an arbitray lattice. 

In this paper, we show that the one-point Green's function is singular in 
the next-to-leading order in the large N expansion. We first evaluate the 
one-point Green's function and the density of state to order 1/N 2 by the 
diagrammatic method for one matrix model. The singular behavior of the 
one-point Green's function in order 1/N 2 has been noticed in the literatures 
[11], but our discussion of this singularity in the context of the two-state 
scattering problem in condensed matter physics may be new. Next we discuss 
the origin of this spurious singularity by the orthogonal polynomial method. 

On the other hand, if we fix iV to be large but finite, and let the energy 
E go to zero, the density of state oscillates and eventually goes to zero. This 
phenomenon is due to energy-level repulsion. The discussion here is some- 
what reminiscent of the double pole encountered in the connected correlation 
function in the Hermitian matrix model when evaluated in the diagrammatic 
approach [2]. In the diagrammatic approach, to calculate the Green's func- 
tions we select diagrams by letting N go to infinity first. We then obtain the 
correlation function by taking the imaginary part of the two-point Green's 
function. The connected correlation function has a double pole as its two 
arguments approach each other. In contrast, in the orthogonal polynomial 
approach, we in effect take the imaginary part of the two-point Green's func- 
tion first, and then let N go to infinity. The connected two-point correlation 
calculated in this way does not have a double pole when its two arguments 
approach each other. However, when we smooth out the short distance oscil- 
lations of the correlation function by averaging it appropriately, we recover 
the double pole obtained in the diagrammatic approach [1]. 

This paper is organized as follows. In section two, the two-state scattering 
problem is formulated as a complex matrix model. We evaluate the one 
point Green's function diagrammatically and discuss the singularity of the 
density of state in d — 0. In section three, we study the complex matrix 
model further using the orthogonal polynomial method. In section four, we 
develop the analysis for general dimensions, using the "lattice of matrices" 
formulation given in [3]. Using the diagrammatic expansion [12,13], we obtain 
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the expression for the one-point Green's function. We compare our results 
with that obtained in previous studies [14-16]. 



2 Matrix model formulation of two-state scat- 
tering 

In the standard Hermitian one-matrix model, the one-point Green's function 
G(z) is defined by 

where the average is taken with the probability distribution P((p) 

P(<p) = exp[-NTrV(<p)] (2.2) 
The N x N matrix tp is Hermitian. In the Gaussian case, we have 

V(<p) = ±Tr<p 2 (2.3) 

More generally, we have, for example, as 

V(<p) = ^Try. 2 + ^TV (2.4) 

For application to disordered systems, the random matrix ip is interpreted as 
the Hamiltonian. 

We start with the simplest case of two-state scattering. The model is 
described by the random Hamiltonian [3], 



taken from the Gaussian distribution P(H), 

p(ff) = l e -iVTr[i(m2if2+m|if 2 2 )+mV+¥'] ^.6) 

z 

The matrices H 1 and H 2 are Hermitian while <p is complex. In [3] this 
Hamiltonian was taken to describe a system with two sectors (C=2). Here 
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we can think of the two sectors as representing as the spin up and spin down 
sectors in a spin-dependent quantum Hall system. We may also think of 
possible applications to the double-layered quantum Hall system [17,18]. 

We now go to the model of (2.6) by letting ml = m\ = oo, so that Hi and 
H 2 are suppressed. This model is for the off-diagonal disorder. This same 
complex matrix model is considered also in [19]. (We will treat the more 
general case with Hi and H 2 non- zero later in this paper.) We have 

h = (°. f) (2-7) 



Notice that there exists a matrix 



I 
-/ 



(2.8) 



such that 

{T,H} = (2.9) 

This implies that if ip is an eigenstate of H with eigenvalue E, then Yip is an 
eigenstate with eigenvalue — E. Thus, eigenvalues of H come in pairs. Due 
to level repulsion, around E — 0, there should be a "hole" of width of order 
1/N in the density of state p(E). As N goes to infinity, this hole disappears 
and p(E) should become smooth. 

To proceed, we calculate the one point Green's function, with P = 
Z~ l exp(—^trH 2 ) where H is given by (2.7). 

G « =< >=< i Tr <i^v» > < 2 ' 10 ' 

To leading order in the large N limit, we obtain easily 

G(z) = Z ~^~ A (2.11) 

using for example the diagrammatic approach of [2]. The density of state 
p(E) = —^lmG(E) is given by the semi-circle law 

p{E) = ^V4^& (2.12) 
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We denote (2.11) by Gq(z) hereafter. 

To leading order, as expected, the density of state p(E) is smooth, without 
any singularity at E — 0. However, if we go to order 1/iV 2 , we will find a 
divergent term in the one point Green's function. Henceforth, we will work 
with the Gaussian distribution. (We expect that the singularity at E — 
will occur also for a non-Gaussian distribution.) 

Using the diagrammatic approach, we can readily evaluate the Green's 
function to order 1/N 2 [2,12,13]. We decompose the self- energy E(z) into 
two parts E a and E&, obtained by breaking the solid line (quark line in the 
terminology of [2]) in the diagrams D a and Dt,, respectively. 

In the simplest two-state scattering case considered here, in which scat- 
tering from 1 to 2 and from 2 to 1 occurs, but not from 1 to 1 or from 2 to 2, 
we see that the numbers n and m appearing in the diagram D a , describing 
the number of rungs in the gluon ladders, must be even. In the diagram D b , 
ni,n,2 and n 3 must be all even, or all odd. We have from diagram D a two 
terms as follows. Denoting the even number of rungs in the ladder by 2n and 
2m, we have 

oo oo 

S a = Yl G An+4m+1 + Yl G* n+4m+1 (2n-l) 

n,m=l n,m=l 

ZU ° (2.13) 



;i - ctr 



The factor of (2n — 1) is due to the number of ways of inserting a cut inside 
the ladder. Similarly, for E 61 , we have 



£ w _ ^ £<4n+4m+4Z+l _^ ^ £,2(2n-l)+2(2m-l)+2(2J-l)+l 
n,m,l=l n,m,l=l 

(1-Gg) s 

For Eb 2 , in which the cut appears inside the ladder, we have 

oo 

S 52 = J2 Gl n+Am+Al+ \2n-l) 

n,m,l=l 
oo 

_l_ (^2(2n-l)+2(2m-l)+2(2«-l)+1^2 n _ 2) 

n,m,l=l 



(2.14) 
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2GKI + GI) 2GI + GI 
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Thus we get 



(1 _ G 4 )4 (1 _ G 4 )3 



2G&1 + Gg) G? 



(2.15) 



(1_ G 4 )4 (i-Gg)s 
Adding the two terms S a and S fe , we have 



(2.16) 



Q7 

s = (i + Gini - giy (2 - 17) 

where Go is the one-point Green's function (2.9) evaluated to leading order in 
the large N limit. After including an extra factor 1/(1 — Gq) for the external 
legs, we obtain the one-point Green's function to order 1/N 2 

Since G (z) goes to (— i) as z — > 0, the factor 1/(1 + Gq) diverges. Using 
1 + Gq = zGq, we have 

G ^ = G ^ + ^2(32,4)5/2 + ( 2 - 19 ) 

Thus, to order 1/N 2 , the Green's function G(z) has a singular imaginary part 
i/(32N 2 z 2 ) as z — > 0. This singularity is related to the reflection symmetry 
(parity symmetry) as mentioned earlier. The density of state diverges like 
p(E) — > 1/ (32/T iV 2 i? 2 ). Apparently, this double pole is too singular, since 
the integral of the density of state should be one. 

We consider next the connected two-point correlation function p 2c (z,w), 
which may be obtained from the connected two-point Green's function 
G2c{z,w) as shown in [1,2]. We have from a diagrammatic analysis [2-4], 

tfG^z.w) = -IA^Log[l - G 2 (z)G 2 (w)] 

' 1 -G(«-' 75 ^-,(^)») (2.20) 



G(z)G(w) w v " G{z)G{wy dz Jy dw 
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This leads to 



1 2^ 2 7/; 2 — A7 2 — Am 2 

N 2 G 2c (z, w) = - 2 £ 2 ™ - 2zw] (2.21) 

A(z 2 - w 2 ) 2 ^Jz 2 - A^/w 2 - A 

The two-point correlation function p2 C (z,w) agrees with the universal 
behavior in the short distance z — > w limit: 

^fc'l-wF? < 2 ' 22 > 

By the equation of motion method, the complex matrix model can be 
studied using a recursion approach order by order in 1/N 2 [11]. Pole terms 
appear in each order of 1/N 2 . The one-point Green's function diverges in 
the k-th order as 

- (2-23) 

Since the connected two-point correlation function in the Hermitian ma- 
trix model does not have a double pole in the orthogonal polynomial ap- 
proach, let us use the orthogonal polynomial approach to investigate whether 
we have the same situation for the one-point Green's function G(z) in the 
limit z — > 0. 



3 Orthogonal polynomial analysis 

Since y? + y? can be regarded as a Hermitian matrix with positive eigenval- 
ues, we can write the joint probability distribution in terms of its positive 
eigenvalues [19,20] 

Pn(£i, ■■■ 1 e N )dsi • • • de N = Ce~ N ^ £l Yi i< j(e i — Sj) 2 dei • ■ ■ de N (3.1) 

with C a normalization constant. The relevant orthogonal polynomial is the 
Laguerre polynomial L n (x), where L (x) = 1, Li(x) = 1—x, and L 2 (x) — 1 — 
2x + x 2 /2. The density of state p(e) and the two-point connected correlation 
function 

P2c(s, £') are given in terms of the kernel K(e,e') [1] as 

p(e)=K(e,e) (3.2) 
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p 2c (e,e') = -[K(e,e')] 2 (3.3) 

where 

i N-l 

K(e,e') = -Y,Me)Me') (3.4) 
iV o 

and 

Me) = e- m ' 2 L n {Ne) (3.5) 

Noting that e = E 2 in our terminology, we have an extra factor of E for the 
density of state since de = 1EdE. At E — 0, all Laguerre polynomials are 
equal to one, so the density of state should vanish at E — due to this extra 
factor E in accordance with the reflection symmetry argument given earlier. 

This result for the density of state at E — seems to contradict the 
result we obtained in the previous section. But there is no contradiction. In 
the previous section, we took the large N limit first, so that the density is 
effectively smoothed over a certain width. Here we take N large but fixed 
and let E go to zero, and obtain p(E) = 0. We see the hole described earlier. 
This non-commutativity of the two limits discussed here is similar to the 
discussion of the two-point connected correlation function of a Hermitian 
matrix model [1, 2]. 

Using the Christoffel-Darboux identity, we have a compact expression for 
the density of state, 

N-l 

P(e) = e-^^LKNe) 

= Ne-x^L^Ne^'n^Ne) - L^Ne^Ne)}. (3.6) 
Since e = E 2 , the density of state p(E) is 

JV-l 

p(E) = 2Ee~ NE2 J2 Ll(NE 2 ) (3.7) 

k=0 

In Fig. 2, this density of state is shown for N = 5 and iV = 10. There 
appears N-oscillations in the density of state, and the first peak near E = 
is finite for iV — > 0. The ratio of the value of the first peak to the second one 
is almost 1.2. In Fig. 2, the dotted line represents the semi-circle behavior in 
the large N limit given by \/A — E 2 / it. When the oscillating part is averaged 
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smoothly, we obtain the correction of the density of state to the semi-circle 
law of order 1/N 2 . 

Ap =< p(E) > --V4 - E 2 (3.8) 

7T 

It may be useful to write the asymptotic expression by the Bessel function. 
Knowing the large N behavior (for E small ) of Laguerre polynomials, we can 
write the preceding in terms of Bessel function. Remarkably, the oscillating 
part near E = in Fig. 2 is approximated by 

p(E) ~ 2NE[J 2 (2NE) + J 2 (2NE)} (3.9) 

By plotting the oscillating curve of this equation, we find that the first peak 
near E = is almost same as the first peak value of (3.7). The ratio of the 
first peak to the second one is also 1.2, which we mentioned before. Now we 
use Hankel's asymptotic expansion of Bessel functions Jo(t), J\(t), for the 
large t = 2NE. 

Mt) = \l^ t [P(0,t)cos(t-^)-Q(0,t)sm(t-^)], 



where 




2 3-7T 3-7T 

Mt) = \l-[P(l,t)cos(t- T )-Q(l,t)sm(t- T )} (3.10) 



P(l,t) = £(-1) 

k=0 



{2t) 2k 



(4/ 2 -l)(4/ 2 -9) 
128t 2 



Q(i,t) = E( 

k=0 



k (J 1 2k + l)_ 
> (2t)2fc+i 



Al 2 -1 
~ 1 

8* 

„ n r(± + Z + m) . . 



m) 



We use the smooth average for the oscillating part by setting < sin (x) > 
=< cos 2 (x) >— \- Then it is easy to find that the density of state becomes 
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The term of order 1/N 2 agrees with the result obtained in the previous 
section by the diagrammatic method. In this approximation, the leading 
term is given simply by - for \/4 — E 2 /ir. Thus we find that the singular 
double pole at order 1/N 2 is in a sense spurious, which appears because we 
took the large N limit first, and it is recovered after the smooth average [1]. 
The density of state does not diverge for E — > 0. Recently the density of 
state of this complex matrix model has been studied in various contexts, and 
related expressions in terms of Bessel functions have been discussed [22]. 

We note that in the case of the Hermitian matrix model, the Bessel func- 
tion is of half-integer order, and consequently we have no poles in the 1/N 
expansion. The double pole 1/x 2 of the connected two-point correlation 
function cancels with a factor sin 2 (x) as x — > [1]. Here we have a different 
situation since sin 2 (a; — |) does not vanish for x — > 0. We must sum the lead- 
ing pole terms of the equation above and it gives eventually a finite result 
for E — > at fixed large N. 



4 Lattice of matrices 

We now extend our analysis to the general (/-dimensional lattice of matrices 
[3]. We place N by iV matrices on the lattice. The gluon propagator is given 
by <r a /3 defined by 

IV1 a/3 

Here M a p, as defined in [3], is a real symmetric matrix whose entries are 
the analogs of m 2 , m 2 ., and m 2 in (2.6). The indices a and f3, which label 
the lattice sites, run from 1 to C, where C denotes the number of sites 
on the lattice. In leading order, the one-point Green's function is given by 
G{z) = jjJ2a9a where g a is determined by 

9a = Fran (42) 

Let us restrict ourselves to a d- dimensional hypercubic lattice, on which 
a quantum particle hops with a hopping matrix is given by a. We introduce 
e(k) by [3] 

cr ay = J2 < a \k > e(&) < k\l > (4.3) 

k 
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For the case of nearest neighbor hopping, we have 



£ W = i + ^£ cosA; « ( 4 - 4 ) 



77T 



Below we will calculate the one-point Green's function for the most general 
case with arbitrary e(k). For specific examples, we often take for simplicity 
the case m 2 = oo and M 2 = 2. 

We will now calculate the one-point Green's function to order 1/N 2 . We 
consider here the general situation defined by some a matrix. In particular, 
for the simple example given in (2.5) we include the diagonal part of the 
Hamiltonian. The self-energy part to order 1/N 2 is obtained from the dia- 
gram of D a and Db, where the number of rungs on the ladders are no longer 
restricted as they were before. We obtain 

S a = G[(- j^^aa? + G[ °_ ^Q2y^Y^Q2 W ( 4 " 5 ) 

Note that repeated indices are not summed unless indicated otherwise. By 
translation invariance this expression is actually independent of a and j3. 
For the self-energy E&, we have two parts E^i and Eb 2 - We obtain for E&i 

nC 2 

Gj:(a^a(a m U(v%aG 2n G 2 ™G 2k = G^^—^Uf (4.6) 

where n, m and k are the number of the gluon propagators. The other part 
E fe2 , obtained by cutting one quark propagator inside the ladder, becomes 

oo 

GEE {o m )p 1 {a n ) ia {a k ) 1 ^)p a G 2n G 2m G 2 ^G 2k 

Pi n,m,j,k 

= G ^ [i T^^ 2) ^ 2{ T^^ 2Ui T^G 2U (4 ' 7) 

Note the rather unusual wasy in which the indices or site labels are arranged. 
Let us check these expression for the case C = 2 discussed earlier. We set 

--(;;) w> 
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The first term in (4.5) becomes 



(_^_) =i(_^ G 2 (4g) 

^l-aG 2 ^ 2 l l-G 2 1 + G 2J 1-G 4 1 ' 



and 



1 — crG 2 a 2 l (l-G 2 ) 2 (1 + G 2 ) 2 ^ (1-G 4 ) 2 K ' 

The factor 1/2 in (4.7) and (4.8) are necessary for a fixed a (which we do not 
sum over.) Thus we get the same result for the self-energy S a = 2G 9 / (1— G 4 ) 3 
as was given earlier (2.11). As for E 6 , we note 

(a n ) a7 = 5« 7 (n = even) 

a7 = a a7 (n = odd) (4-11) 



Using this property, we see that the numbers of gluon propagators n, m, m or 
n + j, m, k should be either all even or all odd. Thus we obtain the previous 
result for S b as given in (4.9) and (4.10). 

We can now immediately go to a rf-dimensional lattice on which matrices 
are placed by inserting the a given in (4.3). Thus, we obtain 

. aG 2 . . aG 2 

aG 2 

= J2 <a \ k><k \~i 7^\ k><k \P> 

k 1-aG 

where we have used < a\k >< k\/3 >= ^e tfc ^~^. Using this diagonalized 
representation, we obtain for the different parts of the self-energy the follow- 
ing expressions: 

1 e k GP 
e k G 
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E« = GE^E^^Hj-z^)] 3 (4-15) 



it 



/3,7 ^ fc 

x ^Ei^^-«)] (4.16) 



where C is the number of lattice points, C = L\ ■ ■ ■ L^. We note that by 
translation invariance £51,62 are in fact independent of a. 

We can easily recover our previous results of course. For the two-site case 
we have k — and k — ir, and f3, 7, = 1, 2 with a fixed a — 1. For instance, 
noting that £a:=o = 1, and £fc =7r = — 1, we get immediately 

G 7 (l + C 6 ) 
(1-G 4 ) 3 



(4.17) 



in agreement with(2.12). Similarly, we find easily that calculated here 
agrees with (2.13). 

We have no divergence for the 3-site lattice, where k = 0, ^ and 
£fc = 1, — |, — |, respectively. Similarly for lattices with odd number of sites. 

In contrast, for a one-dimensional lattice with the number of sites C = 
L = an even integer, we have k = it and e fe=7r = —1, and thus we get a 
divergence when G 2 = —1. 



For L — > 00, we find 



1 ^ cos(£;)G 



2 



V2 



, , cos(A;)G 2J 
= (4,8) 
. r (2G 4 -l + (l-G 4 ) 3 / 2 )(l-y^^) 

J «2 - Lr _ (4.19) 
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where k = 2n(j — 1)/L,j = 1, • • • L For the self-energy parts E&i, Sfe2, the sum 
over (3 and 7 give complications. But the leading singularity of £52 cancels 
exactly the singularity of S a2 in (4.20). The singularity of S bl at G 2 = —1 is 
the same as the singularity of S al . Thus we determine the singularity for one- 
dimensional case at the band center G 2 = — 1 with a single pole divergence 
in order 1/N 2 . The cancellation of the leading singularity of with S a2 
holds for any dimension and coincides with the cancellation found in ref.[5] 
for the lowest Landau level. We have also verified this cancellation by the 
numerical evaluation of the self-energy for large L near G 2 = —1. 

The singularity l/iV 2 |.E| in the density of state as E — * is too strong, 
since the integral of the density of state should be equal to one. In one 
dimension, the off-diagonal disorder case of the tight binding model, in which 
the hopping matrix is real and N = 1 case, is known to have a singularity of 
1 / 1 -£7 1 ( In \E\) 3 in the density of state near E = [14]. In Fig.3, the density 
of state of a finite chain with real hopping matrix for N=l, is evaluated for a 
box distribution. We consider the nearest neighbour random hopping, which 
is represented by the tridiagonal real matrix M, 



M = 





b\ 







61 





h* 








b 2 





h* 








h 






V 



(4.20) 



By the calculation of the density of the eigenvalue of this matrix, we ob- 
tain the curve of the density of state. The random variable 6« is generated 
5000 times, and the histogram of the eigenvalues is evaluated. From these 
calculation, we see the divergent singularity near E = 0. The singularity is 
consistent with 1/E behavior, although it is difficult to see the existence of 
the logarithmic correction from this calculation. 

For the complex hopping case, in which the coupling 6j in (4.20) is a 
complex random number, the density of state shows oscillations similar to 
Fig. 2. We evaluate this complex case in Fig. 4. The first peak near E = is 
finite for L — > 00, where L is the length of the chain. 

It may be also interesting to note that we have the similar behavior in 
different examples. For the sparse random matrix, the density of state shows 
a singularity l/|S|(ln \E\) 2 instead of l/|i?|(ln \E\) 3 [15]. We have also an- 
other example, studied by Brezin, Gross and Itzykson who have obtained the 
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same singularity l/|i?|(ln(|i?|)) 2 in the density of state for the lowest Landau 
level with a random Poisson distribution [16]. 

For the two dimensional case, the self-energy £ ol behaves like [ln(l+G 2 )] 2 
and S a2 like ln(l + G 2 )/(l + G 2 ). This singularity is cancelled with a singu- 
larity in Sfe 2 - We find that the singularity of the self-energy is proportional 
to [ln(l + G 2 )] 2 near the band center. This coincides with the calculation of 
the two-site lowest Landau level result in ref. [5]. 

It is remarkable that we have the same singularity as in the two- state 
(spin up and down) degenerate quantum Hall case studied in [5]. We con- 
sider this a manifestation of the possibility that these two mdoels may same 
universality class. The integral of the density of state with this logarithmic 
square singularity is finite. The exact singularity of the off-diagonal tight 
binding model is not known. Numerical work [9] for the two-state lowest 
Landau level model also shows a singularity at the band center also. More- 
over, the logarithmic square singularity seems to have an applicable region 
near E = according to a recent numerical study for the lowest Landau 
level, although the density of state does not diverge for E — * [23]. The 
state of E = in the two-dimensional case is related to the zero energy 
wave function, for which the Atiyah-Singer index theorem can be applied. In 
this respect, one can perhaps relate the present problem to other interesting 
problems [24,25]. 

For two-dimensional case, we have also evaluated the density of state of 
the off-diagonal disorder (N = 1). We examined the real and complex case 
similar to the one dimensional case. For the real case, the density of state 
seems divergent at E — and consistent with [26]. However, the density 
of state of the complex case, in which nearest neighbour hopping matrix 
element is complex, shows the similar behavior to the complex matrix model 
(Fig. 2). The first peak near E = is finite. 

5 Discussion 

In this paper, we have discussed the singularity of the density of state in 
the large N limit of a complex matrix model as well as its d-dimensional 
lattice generalization. We have found the singularity in order 1/N 2 . We 
have compared our result with known results in d = and d—1. For d — 0, 
the double pole singularity we obtained is reminiscent to the spurious double 
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pole of the two-point correlation function in the Hermitian matrix model. For 
d — 1, we compared our result with the off-diagonal tight binding model. Our 
result is different by the logarithmic factor. For d — 2, we have obtained the 
logarithmic square singularity, which coincides with the result of [5] for the 
A-orbital two-state lowest Landau level quantum Hall model. We interprete 
this coincidence of the singularity as a manifestation of the possibility that 
these models belong to the same universality case. Although our calculation 
is restricted to order 1/A 2 , the result of the singularity of the density of 
state gives a clue of what the true behavior might be. It may be interesting 
to evaluate the Green's functions to order l/N 4 . For the spin degenerate 
two-state quantum Hall system, the numerical similation suggests that there 
are three extended states [9]. In l/N expansion, the singularity (lnE) 2 /N 2 
has been obtained for the density of state [5]. However, as we discussed in 
this paper, this is interpreted as the result after the smooth average. The 
density of state is considered to be finite for E — > 0. It is also important to 
note that the density of state has a nonvanishing value for E — > 0. Then, as 
discussed in [5], the conductivity is exactly given by a xx = e 2 /nh since the 
parameter 9 in [5], defined by 9 = — tan^ 1 (lmG(z)/ReG(z)), becomes —it/2. 
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Figure caption 



Fig. la The diagrams D a and D b of order 1/iV 2 . 

Fig. lb The self-energy S al ,E a2 ,E bl and E 62 , which are obtained from D a 
and D b by cutting the solid lines. 

Fig. 2 The density of state of the complex matrix model for N = 5 (broken 
line) and N = 10 (solid line). The semi-circle law y/A — E 2 /n is also 
shown by a dotted line. 

Fig. 3 The density of state of a finite chain model (L = 10) with off-diagonal 
disorder. The hopping random variable hi is real and obeys the box 
distribution with — 1 < bi < 1. There is a divergence at E — 0. 

Fig. 4 The density of stae of a finite chain model (L = 12) with off-diagonal 
disorder. The hopping random variable hi is complex. The real part 
and the imaginary part both obey the box distribution, which takes 
the value between —1 and 1 uniformly. 
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